The Fermion Determinant, its Modulus and Phase 
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Abstract 



We consider path integration of a fermionic oscillator with a one-parameter family of 



boundary conditions with respect to the time coordinate. The dependence of the fermion 
determinant on these boundary conditions is derived in a closed form with the help of the 
self-adjoint extension of differential operators. The result reveals its crucial dependence on 
them, contrary to the conventional understanding that this dependence becomes negligible 
over sufficiently long time evolution. An example in which such dependence plays a 
significant role is discussed in a model of super symmetric quantum mechanics. 
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X. 

The determinant arising from fermionic path integral over Grassmann variables brings us 
important and critical information on models in field theory. Anomalies fl], [| [| in gauge 
theories are a typical example of this kind. They reveal that the symmetry of the classical 
action breaks down when the fermion fields are integrated out to yield the determinant ||. 
They critically restrict the model in its fermion content Thus any question about the 
determinant deserves careful investigation. The issue of this letter is the dependence of the 
determinant on boundary conditions with respect to the time coordinate. 

It is useful to formulate fermions in field theory by putting the system in a finite spatial 
box. Owing to this treatment the one-particle states can be labeled with an integer. To each 
one-particle state, we assign a fermionic oscillator to distinguish the vacant and occupied states. 
The fermions are then described by a set of the oscillators which interact mutually through their 
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coupling to background bosonic degrees of freedom. For example, in Euclidean path integral 
formalism, a chiral two-component fermion field in a gauge theory is defined by the Lagrangian 



(XT 



where ip a and ip a are anti-commuting fermionic variables which correspond to the creation and 
annihilation operators, respectively, of the a-th one-particle state, v ab is the matrix element of 
the Hamiltonian i[d + iA(x)]a between the a-th and the 6-th one-particle states. Similarly, uj a i, 
arises from the temporal component of the gauge field. To fully comprehend fermions in field 
theory, we will explicitly calculate the determinant of a fermionic oscillator, i.e. a fermion with 
only one one-particle state, with a variety of boundary conditions. 

In the following, we restrict the frequency of the oscillator, z/(r) + iu(r), which is the one- 
component counterpart of the matrices in Eq. ([[]) to be real. This imitates a calculation in 
the temporal gauge (u = 0). Since the frequency is real, two quantum states, vacant and 
occupied, of the oscillator have real and positive-definite transition amplitudes for imaginary 
time evolution. The path integral over anti-periodic configurations is well known to yield the 
sum of the amplitudes of the two states. Thus the determinant is real and positive-definite. 
The periodic boundary condition also yields the sum of their amplitudes, but in which each 
one is weighted by (— ) F , where F denotes the occupation number of the corresponding states 
|| . If the amplitude of the occupied state is larger than that of the vacant, the determinant 
becomes negative. In a previous study [|7| , we considered the determinant with a one-parameter 
family of boundary conditions which includes the periodic and anti-periodic conditions. There, 
our motivation was to clarify the role of a zero-mode in the path integral. We have shown 
the boundary condition which admits it is not the same as the anti-periodic condition and is 
generally different from the periodic one. The determinant becomes zero under this condition. 

Naively, one might think that the dependence of the determinant on boundary conditions 
becomes negligibly small as the time interval goes to infinity. The difference in the values of 
the determinant, however, does not disappear. We extend the calculation to a wider family of 



boundary conditions than we did in the previous study [0] by letting the parameter be complex, 
thus confirming its dependence. We will also provide an example in supersymmetric quantum 
mechanics that shows a crucial dependence on the boundary conditions when analyzed in path 
integration. 

The path integral of the fermionic oscillator is given by 

l-T 



J [dipdip] exp — J dripVip 



(2) 



where T> = d/dr + z/(r) and u(t) is the time-dependent real frequency induced by background 
bosonic degrees of freedom. The integration over the anti-commuting variables, ip and ip, is 
implemented by Berezin [§, || [10, 11], and the result is called the determinant of T>. Since D 
is not a finite-dimensional matrix, its determinant is not simply the product of its eigenvalues. 
Difficulty arises because T> is not hermitian: we use the inner product of vectors, say ip and 4>, 
defined by ((f), ip) = Jq dr(f)*Lp in the Hilbert space composed of square-integrable functions of 
r in the interval [0, T]; the adjoint T>^ = —d/dr + u(r) is then different from T>. We do not 
know a proper setup in the eigenvalue problem of non-hermitian operators that provides us 
with such a useful tool for path integral as a complete orthonormal set of eigenvectors. Even 
the meaning of the eigenvalues of T> is not very clear: under the boundary conditions we deal 
with in the following, the linear manifold on which T> operates is generally different from that 
which results from its operation. 

Fujikawa has proposed a method which is useful in this situation Hl2| . Following him, 
we solve the eigenvalue problems of T>^T> and T>T>^ . Appropriate boundary conditions make 
these operators self-adjoint and non-negative. Then their normalized eigenvectors and (f)^ 



(n = 1,2,3,....) 



have a one-to-one correspondence to each other, 
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This relation holds for all pairs of eigenvectors with a positive eigenvalue, and their spectra, 
the sets of eigenvalues, are identical except a possible difference at zero eigenvalue. Using the 
expansions by these eigenvectors, 

${t) = £a n ^)(r), </S(r) = ^a n ^>(r), (5) 

n n 

we transform the integral measure [dip dip] to [da da] and obtain 

I = N ( Y[[da n da n ]e- ^ = Af [det(P f P)] 1/2 , (6) 

where M is the Jacobian between the two measures, and [det (£>+£>)] 1/2 = n„ V\i is the infinite 
product of the square-root of the eigenvalues. Although the latter is divergent, it gives a finite 
result to I when combined with M. 

The boundary conditions we use are parametrized by one complex variable (3. The conditions 
for T>^T> are given by 

<p(0) + f3(f(T) = 0, (3*V<p(0) + V V (T) = 0, (7) 

and those for T>D^ are given by 

/T0(O) + 0(T) = 0, PV(0)+/3Z>V(r) = 0. (8) 

Note that these conditions keep the one-to-one relation between (p and <p in Eq. (f|). They 
define two linear sub-manifolds in the Hilbert space. These linear manifolds are the domains 
of V^V and W\ respectively^. We can readily verify equations (V^Vipi, </? 2 ) = C^Vi, ^^2) = 
(</?i, T>^T>^p2) for <y9i and ^2 which are arbitrarily chosen in the domain of T^D. Thus X>tX> 
is at least symmetric and non-negative; and so is T)Tft . They are in fact self-adjoint so as 
to have a complete orthonormal basis of eigenvectors, and the expansions in Eq. (|5|) cover all 
possible configurations. The proof is given based on the mathematical theory of the self-adjoint 
extension of symmetric operators |13| ; it will be detailed in a separate publication Jl4| . 



1 We refer to a linear sub-manifold on which a linear operator is defined to act as a domain. 



A different value of (3 provides a different domain for T>^T>. Thus its determinant depends 
on (3 as well as the background v{t), I = Ip[v\. 

We calculate the modulus of I first and its phase next. However, prior to that, we view how 
the determinant has phase, since it is helpful to have a perspective on our calculations. Let us 
start at Eq. (f|). We always take the positive root of y/\ n in this equation. The relative phase 
between cp and <p is then fixed. The consequent determinant [de^T^T*)] 1 / 2 in Eq. (H) becomes 
a product of positive quantities, and it does not induce phase. Hence only the Jacobian Af 
contains the phase of /. 

Now assume temporarily that f3 is real. We can then choose real functions for the eigen- 
vectors and (p^ n \r). These functions can be regarded as the "transformation matrix" 
between iP(t)(iJj(t)) and a n (a n ), according to Eq. fl5|). The Jacobian Af is related to their "de- 
terminant" and is real for this case. At f3 — 1 the determinant is positive-definite, and thus the 
phase in Af is zero independently of v(t). 

For a complex value of /3, the eigenvectors cannot be real and Af becomes complex. To obtain 
its phase, we compare two different measures, one is related to the coefficients in expansion 
(||) where the eigenvectors are solved for the complex (3, and the other is obtained from those 
for (3—1 but with the same frequency v(t). The value of Af differs depending on to which 
measure the original one [dipdip] is transformed. Let ip( n )' and (p^' be normalized eigenvectors 
and [dbdb] be their related measure for the value of (3 of interest; similarly <p^ n \ and [dada] 
are obtained for (3 — 1. These measures are related to each other by [dada] = J\dbdb\ where 

J = [det(y? (n) , ipW)}- 1 x [det(0 (n) *, ^ m ) / *)]" 1 . (9) 

The expressions 

(pto^M') an d (</>(")*, y>M'*) in Eq. (g) are the inner product of the n-th and 
the m-th eigenvectors; each of them is the element of the transformation matrix between the 
two different orthonormal complete bases. The power (—1) in Eq. (Q) comes from the fact that 
the measure is for the integration of Grassmann variables || . Since the bases are complete and 
orthonormal, the transformation matrices between them are unitary, and J is a phase factor. 
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M changes by J when the different bases are used, and thus the phase in M is that of J y which 
is also the phase of /. 

Note that each determinant factor in Eq. (^) is in fact not well-defined separately. The 
eigenvectors have arbitrariness in their phase, but each factor is not invariant under the change 
in the phase. It is the one-to-one correspondence between and 0, Eq. (§), that correlates the 
phase of ip to that of <j) and makes J invariant despite this arbitrariness. 

The modulus of the determinant is readily obtained by slightly extending the calculation 
in the previous study ]7j to the case where the parameter f3 can be complex. To obtain the 
modulus, we define a two-by-two matrix 

M(*\ = ( + PM z ;T) u 2 (z;0) + f3u 2 (z;T) \ 

wipyz) - y p* Vu ^ z . 0) + Vux ( z . T) p*Vu 2 (z; 0) + Vu 2 {z- T) ) ' [W) 

where U\ and u 2 are two linearly independent solutions of a ^-parametrized differential equation 

V^T>Ui(z\ t) = zui(z; r), (11) 

solved with initial conditions Ui(z;0) = 1, Ui(z;0) = dui(z; 0) / dr = 0, u 2 (z;0) = 0, and 
u 2 (z; 0) = 1. The complex parameter z in Eq. (|TTD becomes an eigenvalue of T>^V> if and only 
if some linear combination of u\ and u 2 meets the boundary conditions in Eq. (0). In other 
words, there exists a two-component non-zero vector ji that satisfies [M^(z)]ij / jj = 0. This is 
why we choose the particular form for Mp in Eq. (0). The condition that the z is an eigenvalue 
of V^V is equivalent to detMg(z) = 0. 

From this behavior of Mg, we can prove the identity 

det(X>tx> - z) = ~ ( \n-z \ = det M p (z) 
det(VW -z)~ JJi \\ n -z) det M p {z) 

where the tildes used here distinguish two different operators, matrices and eigenvalues which 
are defined with different frequencies, say u(t) and v{r), but with the same value of the pa- 
rameter (3. Two fractional expressions, far left and far right in ([T2D, are meromorphic functions 
of z. Owing to the behavior of Mp, they have poles and zeros at same values of z, i.e. zeros 



at eigenvalues of Tft T> and poles at those of T^T>. Both functions converge to 1 asymptotically 
as \z\ goes to infinity except along the real positive axis 0]. The ratio of the two meromorphic 
functions is thus an analytic function of z that goes to 1 as z goes to infinity in any direction, 
and is necessarily a constant This concludes the proof of Eq. (|T^). 

The identity (U) means that the ratio |Af[det (£>+£>)] ^(/[det M^O)] 1 / 2 does not depend on 
the z/(t) which is used to define D and Mp. We can thereby obtain the functional dependence 
of the modulus on z/(r) by calculating the determinant of Mp. Putting the solutions at z — 0, 



«i(0;r) 
u 2 (0;r) 



exp 
exp 



- / drV(r') 
Jo 



x h + u(0) J dr'exp 



2 I dr"v{r") 
o 



dr' exp 



2 / dr"v{r") 
o 



(13) 



into the expression of Mg(0) in Eq. ( |T0D and taking its determinant, we obtain 



jV , [(M + (3M 1 )(M + (3*M 1 , 



i 1/2 



(14) 



where Ai = exp (1/2) J T dru(r) , A^i = exp —(1/2) J T drvij) and jV' is a positive constant 
that may depend only on (3. 

Eq. (|i~4|) implies J vanishes for /3 = — A^o 2 (note j\4.i = M.$ l ). This is because the boundary 
conditions with this value of j3 admit a zero mode in the domain. Its effect on the determinant 
calculation has been clarified in [[7]. 

We next proceed to the calculation of the phase, In J . Our plan is to integrate infinitesimal 
variations induced in the phase while the parameter f3 moves continuously from 1 to some value 
of interest. u(r) is kept fixed in this process. We use the formula 51ndet M = TtM~ 1 5M that 
holds for an infinitesimal variation of any matrix M. Applying this formula to the variation 
<51nj7" with J in Eq. (||), and using the fact that a set of eigenvectors is complete, we obtain 

51nJ = ]T [(0^,0^' - (n) ) - (v? (n) ,^ n) ' - V {n) )} = £ [(0 H ,0 ( " ) ') - (<^ (n) y n) ')] , (15) 



2 The same argument is used in Ref. |TT[ to calculate a different determinant. 
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where and 0^ are eigenvectors at (3 and primed ones are those at + 5(3. We put the 
one-to-one correspondence < - n - ) = T>ip( n > / \/Xn into (|T3|) . After integrating by parts and using 
the similar one-to-one relation for 0*™)' and the boundary condition (0), we obtain 

1 



5lnJ = J2 



it 



5 -^^>(T)V^(T)+ 1 -5X n 



(16) 



where 5X n is the variation of the n-th eigenvalue under the infinitesimal change of /3. 5X n is 
calculated by the relation 5A n (^,^W) = {^ n \V^V^') - (V*V<p( n \ (pW) as 



5\ n = -!^(»)*(T)Xty< n >(T) - 5 -^V^ n >(T)^ n \T). 

P P 



(17) 



We then obtain 



^V n >*(T)lV n) (T) - ^p^H*(T)^ (n) (T) 



This expression manifestly shows that 5 hi J is purely imaginary as it should be. 

To sum up the terms in Eq. flT8|), we use the resolvent R z = (T>'T> — z)~ x , which has the 
expression 



E 



i 



^ n \r)ip w *{a) 



(n)*/ 



(19) 



An - Z 

as an integral kernel. We notice readily that contour integration of T>R z /z over z along a 
contour that goes around all eigenvalues of T>^T> clockwise in the complex z-plain gives the 
relevant part of the sum. The contour we use here is one that goes below the real positive axis 
towards the origin from infinity until it passes all the eigenvalues and then goes back moving 
above the real axis. Without changing the value of the integral, we can add to it another 
contour that almost makes a circle at infinity but does not cross the real axis, because the 
value of the integral along the latter contour is zero. We now have the integral along the closed 
contour C (see Figure 1), in which only the pole at the origin contributes to the integral. We 
obtain 

'5(3*, 



* In .7=1 



(3* 



VRo(r,a) 



(c.c.) 



=o-=T 



(20) 



as the result of the sum. 

We have assumed that (3 always stays off the value at which one eigenvalue becomes zero 
and Rq becomes singular; the phase of I is obviously meaningless when / is zero. 

In order to calculate the resolvent Rq in Eq. (PU|), we have to know some mathematical 
details about the self-adjoint extension of differential operators. Following a textbook fl3| , we 



briefly describe a specific extension of the operator TftT> to the extent where our calculation of 



Ro appears convincing; the self-contained description will be given in [14j. The operator that 
is extended to be self-adjoint is D^D defined, however, in a domain D under more restrictive 
boundary conditions than Eq. (|7]). These conditions are <f(0) = <^{T) = 0(0) = <p{T) = 0. T>^T> 
defined in Do is symmetric but is not self-adjoint: the domain of its adjoint is not restricted 
by a boundary condition and is obviously larger than Dq. We have to extend Dq by adding 
two appropriate vectors in order to make T>^T> self-adjoint, that is, the operator itself and its 
adjoint have the same domain. These two vectors have the form 

Wi {z) = U ijUj {z*) - Ui {z) = 1,2), (21) 



where «j is the solution to Eq. (|Tl| ) with an arbitrary but non-real parameter z, and Ui(z*) = 
u*(z) is its complex conjugate. The two-by-two matrix U in Eq. ([H]) has to be "unitary" 
in the sense that the mapping defined by Ui — > UijU* from the two-dimensional vector space 
composed of UiS to that of u*s is unitary. V^V becomes self-adjoint in this extended domain 
D = Do®{ Wl ,w 2 } pj. 

The unitarity of U solely does not determine U, but boundary conditions do. The require- 
ment that any vector in D should meet the boundary conditions in Eq. (0) solves U as 

[U p {z)] t3 = [M p {z%l [M p {z)] ki (22) 

for the parameters (3 and z. The mapping induced by this particular Up(z) is shown to be 
unitary JE4J. The self-adjoint operator V^V that we have used is in fact this extended operator. 



9 



Having described the structure of the domain D, we can now resume the calculation of the 
resolvent. What we need to know in Eq. fl2"U|) is its value only at the end of the time interval. 
Among all the vectors in D, only the two u>jS in Eq. ( PTj ) can contribute to R z at the end. We 
examine the operation of (T>^T> — z) on n>j and find 



(&V - z) Wl (z) = (z* - z)[U p (z)} ijUj (z*), 



(23) 



which enables us to obtain the necessary information on R z . Since zero is not an eigenvalue, 
we can take the limit z — > in Eq. (p3[) to obtain RoUi = Wi, where 



Wj = lim 



Wi(z) 



lim 



{[Up{z)]ijUj{z*) -Ui(z)} 



z^O Z* — Z 2->0 Z* — Z 

and we have used Up(Q) = 1 (see Eq. (p2[)). As an integral kernel, Ro is written as 

R {r,a) = Wi(r)ui(a) + 

i=l,2 



(24) 



(25) 



where UiS are the linear combinations of UiS that satisfy (Ui,Uj) = 5ij, and we have excluded 
terms that do not contribute to the final result. Both Wi and Ui can be explicitly obtained with 
the solutions UiS at z — 0. Plugging these results into (p0| ) and after a long but straightforward 
calculation, we obtain 



(26) 



By integrating this equation under the condition J = 1 at (3 = 1, we finally obtain 

_ r -Mo + /?A4i l 1/2 , 97 , 

The dependence of the phase factor in Eq. (|2"7| ) on the parameter {3 shows the advantage 
of the present calculation where we have allowed (3 to be complex. Let (3 move along the real 
axis, for example, from 1 to —1 and assume the value —M.q 2 is between them. At this point 
the determinant becomes zero, and there will occur an ambiguity in its sign after (3 passing 
the point if we know only its modulus. This ambiguity is solved since we can use a path that 
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circumvents the point along an circle with infinitesimal radius e as (3 = — A4q + ee td where 9 
moves from to it. The determinant becomes negative after (3 passes the value. 
Combining the modulus (|H|) and the phase factor (P7|), we obtain 

Ip[v] =W{M Q + (3Mi). (28) 

M.q and M.\ are in fact the transition amplitudes for the vacant and the occupied states, 
respectively. The result in Eq. (p8|) is not only consistent with the known case of the anti- 
periodic {p = 1) and the periodic (/3 = —1) boundary conditions, but also applies to any 
complex value of f3. Af' is an overall normalization and will not affect the following discussion. 

Although we have considered only a single fermionic oscillator, our result indicates that 
the similar dependence is expected to appear in fermion determinants in field theory, perhaps 
in a more complicated manner. At least for a simple case in which there are a finite number 
of fermionic oscillators that share a common value of (3 for their boundary conditions, we 
can safely conclude, without seeing anomalies, that the determinant is the trace of the time 
evolution operator in which the amplitude of each state is weighted by j3 F , according to its 
fermion number F. It is hard to imagine that the difference of the determinant caused by such 
a non-trivial weighting disappears even in the limit T — > oo. Thus we have to be careful which 
boundary condition we use in the calculations. 

In order to see how the boundary conditions affect the result, let us examine an illuminating 



example — the supersymmetric double- well quantum mechanics fTS ]. With q and p, the coor- 
dinate and momentum of bosonic degrees of freedom, ip^ and ip, the creation and annihilation 
operators for fermion, we write the Hamiltonian of this model as 

H = \{Qtf) = \{V 2 + W{ q f) + (29) 

where Q = (p + iW(q))ip, = (p — iW(q))ip^ and W is chosen to be W(q) = q{l — gq) with 
a coupling constant g. What will happen if one analyzes the model by path integral without 
paying attention to the boundary conditions? One may try to integrate the fermion first in a 
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given bosonic background and then to integrate the bosonic part weighting the bosonic measure 
with the fermion determinant, which is given by Eq. (Q) with v(t) = 1 — 2gq{r). Under the 
periodic boundary condition, the determinant changes the sign when the background varies 
between q(r) and 1/g — q(r), for which the bosonic potential W(q) 2 takes the same value. The 
contribution of the two configurations cancels each other out, and the path integral becomes 
zero. One may conclude that the model is ill-defined by interpreting this result as the statistical 
trace being zero. 

We know, however, that the model is well-defined, and further it is an instructive example of 



the dynamical supersymmetry breaking ||15|| . The path integral over the periodic configurations 
is not the statistical trace, but the regularized Witten index Tr(— ) F e~ TH [fL6[ , |(J. The fact that 
it is zero merely implies the possibility that two perturbative zero-energy states, one localized 
at q = with F = (F = ip'i/i) and the other localized at q — 1/g with F — 1, may be lifted 
in pairs by non-perturbative effect of the interaction, and the supersymmetry may be broken. 
It indeed happens in this model by the valley-instanton effect [17] . 



This example shows that in any model which exhibits dynamical supersymmetry breaking, 
its fermion determinant needs to have non-trivial dependence on the boundary conditions. The 
path integral has to vanish in the case of the periodic fermion fields so that the Witten index 
is zero, while it should provide a non-vanishing value for anti-periodic fermion fields in order 
that the partition function is not zero. Thus the careful study of the dependence of fermion 
determinants on the boundary conditions is important in the search for a model in which 
supersymmetry breaks dynamically. 
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Figure caption: Figure 1. The schematic drawing of the contour C in the complex z-plane 
of VR z /z. C is a closed contour made of a contour that goes around all eigenvalues of V^V 
clockwise and another that goes around at infinity counterclockwise, as explained in the text. 
The crosses indicate positions of the poles of VR z /z. 
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